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Reaction-diffusion model for pattern formation in E.coli swarming colonies with slime. 
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A new experimental colonial pattern and pattern transition observed in E.coli MG1655 swarming 
cells grown on semi-solid agar are described. We present a reaction-diffusion model that, taking 
into account the slime generated by these cells and its influence on the bacterial differentiation and 
motion, reproduces the pattern and successfully predicts the observed changes when the colonial col- 
lective motility is limited. In spite of having small non-hyperflagellated swarming cells, under these 
experimental conditions E. coli MG1655 can very rapidly colonize a surface, with a low branching 
rate, thanks to a strong fluid production and a local incremented density of motile, lubricating cells. 

PACS numbers: 87.18.Hf 87.18.Bb 87.18.Ed 87.80.Vt 



I. INTRODUCTION 

The macroscopic pattern exhibited by a bacterial 
colony and its dependence on certain environmental pa- 
rameters can give us some clues about the coordinated 
colonization strategy followed by the community of cells. 
The biological interest in this interdisciplinary area is 
in pointing out under controlled laboratory conditions 
the cooperative mechanisms (intercellular interactions, 
motility and communications) that these bacteria, which 
are traditionally considered as solitary life forms, may 
have developed to adapt to changing environments. It 
is an extremely fertile and interesting area for close col- 
laboration between physicists and biologists because it 
helps both to understand the transition from individual 
(uni-cellular, a bacteria) to collective (multi-cellular, a 
colony) behavior pj. 

Bacterial colony pattern formation on semi-solid agar 
surfaces has been studied extensively by microbiologists 
and physicists QSHHSSHIl. 

Through these works 
undertaken with colonies arising from different species 
and strains of bacteria, the following common features 
have been described: Bacteria can swim in liquid medium 
without difficulty but in environments with adverse con- 
ditions for the swimming motility, they need to develop 
mechanisms to become more motile. On semi-solid agar 
medium, where the viscosity is high and the motility of 
short swimmer cells very low, some may differentiate into 
very elongated, multinucleated andprofusely flagellated 
swarm cells that can move easier p|. The initiation of 
swarm differentiation seems to be strictly correlated with 
physico-chemical factors such as surface contact and quo- 
rum sensing response (cell density sensing mechanism) 
P| . Swarming cells have the ability to extract water from 
the agar and produce a lubrication fluid (slime). The 
flagellum driven motility of swarmer cells together with 
the extracellular slime helps to overcome the surface fric- 
tion After some time migrating swarmer cells have 
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been observed to cease movement, septate and produce 
groups of swimmer cells (de-differentiation process) 
The chemotaxis system is essential for swarming motil- 
ity in E. coli 0. With the chemotaxis mechanism bac- 
teria can orient their motion in response to a gradient 
(positive or negative) of a certain chemical field (a nu- 
trient or a field produced by the bacterial cells such as 
chemical signals or a pH change) 0. Based on these 
experimental observations great effort has been devoted 
to model the cooperative behavior of a cell community 

mm mm p. 

Here we describe the new pattern and pattern transi- 
tions observed in a swarming colony of E. coli MG1655 
a laboratory domesticated wild type strain that 
produces a lubricating fluid or slime: when the colo- 
nial motility is limited this colony expands irregularly 
with few relatively thick, dense, branches and very rare 
secondary branches, whereas for conditions of improved 
colonial motility the pattern is compact and round with 
structures of higher cell densities. To study the coordi- 
nated, self-organized, colonization scheme of this colony 
and based on experimental microscopic observations, we 
will define a mathematical model of the cooperative be- 
havior of this colony which is able to predict the pattern 
transitions when certain control parameters are changed. 

This work is organized as follows: in Section [H] we 
will introduce the experimental results and their impli- 
cations for the non-expert reader. Then in Section IIIII 
our hypothesis for the collective dynamics will be ex- 
plained, followed by a mathematical formulation. Results 
of the computer simulation in comparisons with exper- 
iments will be posed and discussed in Section ITVl The 
summary and conclusions are given in Section Ivl 



II. EXPERIMENTAL OBSERVATIONS 

The development of an E. coli MG1655 colony on a 
semi-solid surface (0.5% "Difco" Agar) under certain nu- 
tritional (1% "Difco" Tryptone, 0.5% "Difco" yeast ex- 
tract, 0.5% NaCl and 0.5% D-(+)-glucose) and environ- 
mental (at 37°C and on average 22% relative humidity) 
conditions has been studied in detail with round and 
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square laboratory plates of 9 and 12 cm respectively. 
We inoculated a 2.5 fil liquid drop of stationary phase 
bacteria cultivated overnight (in LB, 0], at 37°C). The 
following experimental observations were made: 

• The colony expands on the surface in a compact 
and round shape with few (typically 6 to 8) wide 
branches. The colonization of the plate takes place 
very rapidly: at 37° C a plate of 9 cm diameter is 
colonized in roughly one day, see Fig. ^a) . 

• It was observed that the fronts of two colonies in- 
oculated on the same plate never merge unless they 
meet tangentially. Branches do not simply grow ra- 
dially from their inoculation point, but avoid this 
intercolonial demarcation region, see Fig. ^b) . 

• The term "swarming" is used in the literature to de- 
scribe the active,flagellar-dependent surface motil- 
ity of bacteria y|. To confirm the relevance of 
flagella in the mechanisms of surface translocation 
of this colony we have constructed by PI trans- 
duction Ha a flagellar-defective mutant, MG1655 
BhDr.Kan 17]. The BhDC operon encodes a key 
master regulator of the hierarchical system control- 
ling the synthesis of the bacterial flagellum in E. 
coli 0| . Having this operon inactivated the strain 
MG1655 iJiiD::Kan can thus not synthesize the bac- 
terial flagella. The wild type and the flagellar- 
defective mutant were inoculated simultaneously. 
After one day the mutant had initiated expansion 
over the surface, due to the combined influence of 
the expansive forces of bacterial growth (popula- 
tion pressure) and the fluid production. But, in 
contrast with the wild type MG1655, it failed to 
continue the rapid surface colonization, see Fig. 
Hfc). The slow passive form of surface transloca- 
tion displayed by the mutant is called "spreading" 
|19|. The same pattern was also observed when 
the flagellar-defective mutant was inoculated alone 
on a plate and for lower agar concentration (data 
not shown). Notice also the existence of a clear 
asymmetric intercolonial demarcation region which 
has been avoided by the fast expanding wild type 
colony MG1655. This experiment shows thus that 
flagellar motility is essential for the surface motil- 
ity of MG1655 and that we are dealing with an ac- 
tive process of flagellar-dependent swarming motil- 
ity that is able to avoid a region. 

• The macroscopic characteristics of the experimen- 
tal patterns are extremely sensitive to humidity 
which in turn limits the fluid production. When 
the fluid production is high the pattern is almost 
compact with circular envelope and has few wide 
branches with cells deposits and channels struc- 
tures within. Whereas when the fluid production is 
limited the colony expands irregularly, with long, 
relatively thick branches with high cell densities 



(see Fig. ^d)), very few secondary branches and 
greater gaps (see Fig. EJb)(d) and (f)). 

• There is a transparent surrounding envelope of 
slime produced by the colony, see the lower arrow 
in Fig. I^a), and channel- like structures within the 
colony surface, see the upper arrow in Fig. [2Ja). 

• Along the edge in contact with the slime, one can 
observe a thin layer of quiescent cells apparently 
packed in an ordered way. See the lighter line of 
packed swarmer cells in Fig. |2Jb), pushed by the 
turbulent collective motion of groups of swarmer 
cells in the interior region. 

• In this swarming pattern we found, both in the cen- 
ter and the periphery of the pattern, mononucle- 
ated active cells, see Fig.EJa), from 3 to 6 /urn and 
also with 4 to 6 flagella, see Fig.[3{b). The swim- 
mer cells of E. coli MG1655 found in other patterns 
with interstitial submerging are roughly 2 /im long 
with 4 to 6 flagella (data not shown). Therefore, 
under these experimental conditions the swarmer 
cells of E. coli MG1655 do not seem to undergo 
a significant change in their cellular characteristics 
for surface translocation. 

In order to choose an appropriate model certain fea- 
tures must be emphasized from the experimental observa- 
tions: In spite of the small, non-hyperflagellated swarm- 
ing cells, the colony expands very rapidly on the agar 
showing a strong sensitivity to the surface humidity con- 
ditions, we therefore conclude that the lubricating fluid 
is relevant for the successful expansion of the colony. 
Furthermore, we believe that the channel-like three di- 
mensional structures are accumulations of bacteria in re- 
gions with higher concentration of fluid. As for the inter- 
colonial demarcation region, our hypothesis is that the 
motion of motile cells is driven by negative chemotaxis 
avoiding the regions with high concentrations of cells and 
the substances produced by them. Notice in particu- 
lar the wide intercolonial demarcation zone in Fig. ^c). 
Thus despite the fact that the mutant colony continues to 
grow (notice the long bacterial filament) which indicates 
the presence of nutrient in this zone, the swarming wild 
type colony avoids this region. This suggests that the 
demarcation zone cannot solely be due to nutrient deple- 
tion and there must be some chemotaxis in response to 
chemical signaling, nutrient level, pH variations, etc. 

There are a few other reported examples of colonial 
pattern formation in the literature, most of them im- 
ply the splitting formation of tips and branches but lead 
to different patterns since the mechanisms of advance 
and self-organization are also different. We propose to 
describe this dynamical system by means of a reaction- 
diffusion system, a standard tool in the field of bacte- 
rial colony pattern formation, including well known and 
widely used features such as lubrication and chemotaxis. 
The collective diffusion term of the colony must account 
for the new self-organization mechanisms observed in our 



3 



(a) (b) 




FIG. 1: (a) (Color online) Surface colonial pattern ob- 
served in E. coli MG1655 strain in semi-solid agar. (b) 
(Color online) Two colonies grown on the same plate avoid 
each other, (c) (Color online) Comparison of the swarm- 
ing flagellar-dependent surface colonization of the wild type 
MG1655 inoculated on the left, with the spreading growth of 
the flagella-defective mutant MG1655 flhD::Kan inoculated 
on the right. The fast advancing swarming colony avoids the 
slowly growing one leaving an intercolonial demarcation re- 
gion, (d) (Color online) Close view (box of 2 cm side) of a ram- 
ified pattern when the experimental conditions are changed 
(grown at lower temperature 30° C) showing relatively thick 
branches with increased cell densities. 



(a) (b) 




FIG. 2: (a) (Color online) Close view of colony edge with 
surrounding slime (lower arrow) and channel-like structure 
(upper arrow). (Bar=5 mm), (b) (Color online) Swarming 
cells at the edge of the colony, bounded with a thin layer of 
quiescent and packed swarming cells. (Bar=50 /im). 
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FIG. 3: (a) (Color online) Mononucleated E. coli MG1655 
swarmer cells, stained by DAPI fluorescence dye and visual- 
ized by phase contrast microscopy. (Bar=5 /im). (b) (Color 
online) Flagella of swarming cells visualized with the Saizawa 
and Sugawara procedure Q. (Bar=5 /im). 

experiments: we have seen that an increased density of 
cells and slime in the tips of the colony helps the colony 
to advance, being less sensitive to local irregularities and 
showing thus, in conditions of adverse motility, a poorly 
branched pattern. 



III. REACTION-DIFFUSION MODELING 

Previous studies on pattern formation of other lubri- 
cating bacteria such as P. dendritiformis var. dendron 
have confirmed that reaction-diffusion models with a lu- 
bricating fluid and chemotactic motion can successfully 
describe the development of branching structures 
However, when the colonial motility is decreased these 
models predict a pattern transition from a compact to 
a highly ramified pattern with thin branches, which is 
different from the one observed here. 

We will next construct a mathematical model of the 
reaction-diffusion type, to describe our experimental ob- 
servations. The cells in the colony can be in a normal 
static phase S or in a differentiated motile state V where 
the cells produce liquid and use their flagella to move. 
Let U be the nutrient consumed by cells in their repro- 
duction, see Eq. ©. We include the fluid substance C 
that is excreted at a rate (3 by motile cells and disap- 
pears with time at a rate /i c (being reabsorbed into the 
agar medium or consumed or degraded by the cells), see 
Eq. The initial condition is the circular inoculum 

(with radius R) of a drop of liquid C with bacteria S in 
the static state in a plate with nutrient Uq. When the 
concentration of fluid is sufficiently high, greater than a 
threshold C c , these static cells reproduce at a rate A into 
more static cells (in our model their reproduction rate is 
increased in liquid media), see Eq. (J2J. When the fluid 
level is too low, i.e. drops below the critical value C c , 
and bacteria are in contact with the surface the swarm- 
ing behavior is triggered: the new cells are motile cells 
V that are able to produce fluid and move, see Eq. QJ. 
With time, motile cells return to the static state at a rate 
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/j. The diffusion term in Eq. can be expressed as VF 
with F = Vvy a flux of motile cells of density V and 
with velocity vy proportional to V, C and the gradient 
of V. This velocity represents the fact that motile cells 
move faster for higher concentrations of motile cells V 
and lubrication fluid C, and that the direction of their 
motion is actively defined by negative chemotaxis moving 
towards regions of lower motile cell densities V and sub- 
stances produced by them (this negative gradient term 
could alternatively describe a population pressure in the 
case of passive translocation). Let D v be the diffusion 
parameter of motile cells. The nutrient and fluid fields 
also diffuse with diffusion constants D u and D c respec- 
tively. To introduce the local irregularities of the agar we 
include a local random viscosity term rj(x, y) — £ (with £ 
a random number of uniform distribution in the interval 
[0, 1]) which is constant in time and affects the nutrient 
diffusion. Our adimcnsional model reads: 



dV 

— = +XUSCQ (Cc 



-c) 



liV + D v V ■ [V 2 CVV] , (1) 



dS_ 
~dt 



= +XUSCe {C -c c ) + fJ>V, 



dU 
~dt 



= -XUSC + D uV (x, y)AU, 



dc 

dt 



+(3V - n c C + D C AC, 



(2) 



(3) 



(4) 



where O is the Heaviside step function (other thresh- 
olding functions lead to similar results). On swarming 
experiments a lag phase is always observed after inoc- 
ulation when bacteria proliferate by cell division at the 
central spot without any migration. Once the colony 
has reached some threshold, the swarming process starts, 
i.e. rapid surface migration is preceded by a cell density- 
dependent lag period The thresholding condition of 
our model allows us to describe the biological colonial 
switch between the swimming and swarming state and 
this aforementioned initial time lag between inoculation 
and colony spreading. 

We have solved this system numerically using a fi- 
nite difference scheme and an alternating direction tech- 
nique for the diffusion term. We integrate inside a cir- 
cular region of radius 112.5 and set the initial condi- 
tions to U(x,y;t = 0) = Uo, V(x,y;t = 0) = 0, 
C(x, y;t = 0) = £(x, y)Q[ R 2_ {x 2 +y 2 )] and S(x, y;t = 0) = 
£(x, y)0[^2_( :E 2 +J( 2)] with R = 5 and £ a random num- 
ber of uniform distribution in the interval [0, 1]. We fix 
fi = 0.15, /a c = 0.5, D c = 0.015, D u = D v = 1.5 and 
C c = 1.5 and let other parameters vary depending on 
the experiment. 

If we set X = 1, Uo = 2 and [3 = 1.3 we reproduce 
the experimental pattern. Figures 0] (a) and (c) display 
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FIG. 4: Numerical simulations results of reaction-diffusion 
growth with slime, when A = 1, Uo = 2 and /3 — 1.3. Local 
concentration of motile and static cells V + S (a) pattern from 
single inoculation, shown at t — 650, (b) same study case, 
local density of fluid C, (c) pattern from double inoculation, 
shown at t = 465. 



the local bacterial density V + S at a certain time for one 
and two inoculi respectively. The color coding ranges be- 
tween white for zero density and black for values equal to 
or higher than 5. The inner black structures correspond 
to regions where temporary high concentrations of fluid 
led to higher reproduction rates and an increment in the 
static bacterial density. The demarcation region induced 
by the gradient of the diffusion term can also be distin- 
guished. For the study case (a) , we show the local density 
of fluid C at the same given time, see Fig. 0] (b). This 
fluid will disappear with time (it is reabsorbed into the 
agar or used by the bacteria behind the front). 



IV. NUMERICAL RESULTS AND 
EXPERIMENTAL VALIDATION 



The model predicts a pattern transition from compact 
to branched growth as the colony motility is limitted 
due to a reduction in any of the factors in the collec- 
tive diffusion term in Eq. Q: D v , V or C. For instance 
if the reproduction rate (A), the nutrient content (Uo) 
or the possibility to produce slime using water from the 
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FIG. 5: Numerical (left) and experimental (right) results 
when the production of fluid is limited, (a) A = 0.6, Uo = 2, 
j3 = 1.3 at t = 875 (b) (Color online) colony of size 6x6 cm 2 
approximately, grown on a square plate with 0.5% glucose at 
lower temperature (30°C) (c) A = 1, Uo = 1.5, f3 = 1.3 at 
t = 600 (d) (Color online) grown with 0.3% glucose at 37°C 
(e) A = 1, Uo = 1.5, f3 = 1 at t = 750 (f) (Color online) grown 
with 0.3% glucose at 37° C, when the surface lost 1 g of water 
due to direct heating. 



agar (/?) are decreased separately the pattern becomes 
less compact and with an irregular envelope, with rela- 
tively thick branches and very rare secondary branches. 
These numerical results are shown in Fig.[5fa),(c),(e) and 
compared with experimental results where the collective 
motility is reduced by lowering the temperature (b), the 
glucose content (d) or the surface humidity (f) respec- 
tively. Experimentally the most critical parameter was 
the humidity: variations of only 5% of the surface humid- 
ity had a significant impact on the emerging patterns. 

Although the simple pattern analogy is not sufficient to 
conclude that the model is correct, predicting the behav- 
ior when several control parameters are changed is a good 
tool to validate it. Under the circumstances reported 
here the parameter that had the strongest impact on the 
emerging pattern was indeed the humidity. But under 



different experimental conditions other factors, such as 
local cell densities and quorum sensing or type and quan- 
tity of nutrient supply required for the flagellar synthesis 
and activation, may be considered as essential. 



V. CONCLUSIONS 

Branching is due to the sensitivity of the system to 
local irregularities whose small perturbing effect is en- 
larged due to the non-linearity of the diffusion term de- 
scribing the cooperative spreading of the colony. Branch- 
ing patterns have also been seen in other bacteria such 
as B. subtilis, 20], and in P. dendritiformis var. den- 
dron (genera Paenibacillus of B. subtilis), [Uj]- In these 
colonies, when the motility of the colony is strongly lim- 
ited, the advancing front is extremely sensitive to local 
irregularities and shows very frequent tip splitting lead- 
ing to very ramified, thin branches. In some cases (like 
the so-called Dense Branching Morphology) the distri- 
bution of branch lengths is even exponential, suggesting 
that the tip-splitting of branches takes place at random 

In our study case when the motility is reduced, the 
density of motile lubricating cells in the advancing front 
is incremented, the fluid level increases also and the ad- 
vancing branches arc thick, with low branching rate. This 
mechanism of response of the colony to adverse condi- 
tions has been implemented mathematically as follows: 
a thresholding condition on the fluid level launches the 
production of lubricating swarming cells V and the re- 
production rate is increased in fluid rich regions. Since 
the collective colonial motility expressed by the diffusion 
term shows a strong dependence on the density of motile 
cells V, the very dense branches are less sensitive to local 
irregularities and show thus a lower branching rate. 

Studies such as this one are aimed at understanding 
the principles underlying self-organized pattern forma- 
tion processes and establishing a systematic correspon- 
dence between the observed macroscopic colonial pat- 
terns and the putative microscopic collective physics. 
The new biological pattern has indeed given us some 
hints about the underlying collective physical principles 
or colonial strategy. 

This strategy, or part of it, may be common to other 
bacterial patterns. For instance branching patterns have 
been reported recently on a colony of B. subtilis where 
potassium ions were shown to increase the production of 
lubricating fluid with surfactin and thus increment the 
motility substantially without formation of flagella (not 
swarming)!^. When the motility is reduced due to a 
low fluid production the colony shows also thick branches 
with relatively low branching rate. A fluid preceding each 
dendritic finger of growth can be seen. From analysis 
of this information alone, we conjecture that a similar 
mechanism may take place in this case: under adverse 
conditions the production of lubricating cells would be in- 
creased, reproduction may also be more effective in fluid 
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rich regions and in turn these branches would have higher 
density showing again a low sensitivity to local irregular- 
ities and low branching rate. The population pressure 
would also induce here a negative gradient term. 

Swarming in E. coli has been previously reported on 
other strains that swarm due to a differentiation from the 
swimming state to a swarming-adapted strongly hyper- 
flagellated, multinucleated state |2l|. This present work 
reports on the swarming ability of another E. coli strain, 
the wild type MG1655, that shows a greater variety of 
patterns. Instead of cell hyperflagellation, the strategy 
chosen for colonization is a strong fluid production and 
incremented density of motile cells, the result is a fast 
advancing colony with low branching rate. Since E. coli 
is a genetically well characterized organism, the exper- 
iments reported here can be used for further studies of 
the genetic expression associated to the collective, self- 



organized, swarming pattern formation process. 
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